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Abstract. In this paper, a weak formulation of the discontinuous vari- 
able coefficient Poisson equation with interfacial jumps is studied. The 
existence, uniqueness and regularity of solutions of this problem are ob- 
tained. It is shown that the application of the Ghost Fluid Method 
by Fedkiw, Kang, and Liu to this problem in can be obtained in 
a natural way through discretization of the weak formulation. An ab- 
stract framework is given for proving the convergence of finite difference 
methods derived from a weak problem, and as a consequence, the Ghost 
Fluid Method is proven to be convergent. 



1. Introduction 

This paper proves the convergence of the finite difference method intro- 
duced in for the Poisson equation with discontinuous coefficients and 
given interfacial jumps. Based on the Ghost Fluid Method (GFM) of [Q], 
see also [^], the finite difference method is simple, efficient and robust. 
One of the novelties, and advantages, of the method is the arm-by-arm split- 
ting technique which makes the method in multi-dimension as simple as in 
one-dimension. Another advantage is that the resulting linear system of the 
method is the same as the linear system obtained from the simplest standard 
five point stencil finite difference method for the Poisson equation without 
discontinuous coefficients or given interfacial jumps. Therefore the result- 
ing linear system is symmetric and positive definite and can be efficiently 
solved. Inherited from the GFM, this finite difference method captures the 
sharp solution profile at the interfaces without smearing. This is a necessity 
in the development of effective schemes for problems involving interfaces. A 
variety of other approaches to problems with interfaces have been taken, , 
Q, [10 1, |11|, |jl2|. An important application of the method here is for 



Hele-Shaw flow, see |6[. 

The convergence proof starts with the formulation of the problem, in 
terms of a uniformly elliptic bilinear form. Existence and uniqueness of 
a weak solution v follow immediately using elementary functional analysis. 
The solution space is the same as Hq except with a different, but equivalent, 
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inner product induced by the bilinear form. Discretizing this weak formu- 
lation in a natural way results a discrete weak problem, which is equivalent 
to the finite different method in As expected, the discrete weak formu- 
lation inherits the structural conditions (uniform boundedness, extension, 
consistency) in the discrete sense, hence existence, uniqueness, and uniform 
boundedness of the family of discrete solution hold. 

We then provide an abstract framework for the convergence proof. In 
addition to the uniform structural conditions for the weak problem and its 
discretizations, we postulate the existence of a uniformly bounded family of 
extension operators from the discrete spaces to the solution space, which 
satisfy a strong approximation property. This natural assumption implies 
that as the grid size h goes to zero the image of the discrete space fills out the 
entire solution space. Finally, we impose weak consistency of the discrete 
problem with the original problem. Under these assumptions, the sequence 
of extended weak solutions T^{v^) converges weakly to the solution v in Hq. 
Guided by the abstract framework, we establish the converge of the finite 
difference scheme of |9|. 

There is a similarity between the standard convergence proof of finite el- 
ement methods and our approach. Both of them use structural conditions. 
For finite element methods, the structural conditions and weak consistency 
are inherited directly from the weak problem for the PDF because the dis- 
crete bilinear forms are obtained by restriction to finite dimensional sub- 
spaces. Cea's Lemma then says that the extended discrete solution is the 
closest function in the finite dimensional subspace to the true solution. This 
reduces convergence and error estimation to a problem in approximation 
theory. For finite difference methods, further approximations are made so 
that the discrete problem can not be obtained by restriction to a finite di- 
mensional subspace. In our case, the difference between the finite difference 
scheme and the finite element method is that point values of the coefficients 
are used rather than cell averages. Because the coefficients in our prob- 
lem are discontinuous, the solution lies only in Hq and not fl H^, and 
therefore, we obtain convergence, but without a rate. 

2. Equations and Weak Formulation 

Consider a bounded domain, C M"", with smooth boundary, di}. Let 
r C be a smooth interface of co-dimension n — 1, represented by the zero 
level-set of a smooth function (p{x), which is a signed distance function of 
the interface locally. We assume that (j) divides into disjoint subdomains, 
il^ = {(f) < 0} and O"^ = {(p > 0}, with dU^ = T. Thus, we may write 

= il+ U U P. The unit normal vector of the interface is n = V<^/|V</>|, 
for ^(x) = 0, pointing from 0" to fi"*". 

We seek solutions of the variable coefficient Poisson equation away from 
the interface given by 



(2.1a) 



V- (/3(x)Vn(x)) = /(x), xen\r, 
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in which x = (xi, . . . ,Xn) denotes the spatial variables and V is the gra- 
dient operator. The coefficient /3(x) is assumed to be a positive definite, 
symmetric nxn matrix, the components of which are continuously differen- 
tiable on the closure of each disjoint subdomain, 0,^ and fi^, but they may 
be discontinuous across the interface T. It follows that there are positive 
constants m < M with m Id < /3(x) < M Id, where Id stands for the nxn 
identity matrix. We suppose that on the interface, f3 assumes the limiting 
values from within O". The right-hand side /(x) is assumed to lie in I^iyt). 

Given functions a and h along the interface F, we prescribe the jump 
conditions 

J [n]p(x) =u+(x) -n-(x) =a(x), 
^'■'^^ \ [{^uUr (x) ^ (/3n)+(x) - {(3u)~{x) = b{x), ^ ^ ' ' 

Note that {Pu)n = n ■ P'Vu, and the "it" subscripts refer to limits taken 
from within the subdomains O^. 

Finally, we prescribe boundary conditions 

(2.1c) n(x) = g{x), x E dft, 

for a given function g on the boundary. 

We are going to use the usual Sobolev spaces IIq{CI) and II^{Q,), but 
instead of the usual inner product we choose one which is better suited to 
our problem: 

(2.2) B[u,v] = / pVu-Vv. 



This induces a norm on Hq{Q,) which is equivalent to the usual one, thanks 
to the Poincare inequality and the uniform bounds for the coefficient matrix. 

Let i?r and Rqq denote the restriction operators from 11^(0.) to L^(F) 
and L^((9il), respectively. Throughout this section, we shall always assume 
that our boundary data a, b are the restrictions of functions a, b € IIq{Q), 
and that g is the restriction of a function g G H^{Q,): 

(2.3) a = i?r(S), b = Rr(b), and g = Ranig)- 

To simplify the notation, from now on we will drop the tildes. 

First let us consider the version of problem of (2.1a), ( p.lbj ), (2.1c) with 
homogeneous boundary conditions for u: 



(2.4) 



V • (/?(x)Vu(x)) = /(x), x€0\F 

[n]p (x) = 0, X € F 

[(/3u)„]r(x) = 6(x), xGF 

u{x) = 0, X G dQ. 
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Definition 2.1. A function u G Hq{Q) is a weak solution of ( |2.4D if it 

satisfies 

(2.5) -B[u, i;]= J fij+ Jv- (b^n) , 

n n- 

for all V G -f^o(^)- 

A classical solution of (2.4), u\q± € (7^(0='=), is easily seen to satisfy (2.5). 
The boundary condition [u]r = is understood with the aid of i?r, and it 
is fulfilled since u G Hq{Q). 

Theorem 2.1. Given f € L'^{Q) and b € Hq{^), there exists a unique weak 
solution of (|2.5| ) in Hq{^). 

Proof. The right-hand side of (2.5) 



n 



is a continuous linear functional on Hq{^1,). By the Riesz representation 
theorem, there exists a unique u G -ffo(f^) such that —B[u,ip] = F{ip), for 
all V G -?^o(^)- O 

Next we reduce the general case (|2T^ ), (|2lb| ), ( ^1^ to the homo geneous 
case. We will construct a unique solution of the problem in the class 

H{a,g) = {u:u- g + axn- G -f^o(^)}, 

in which Xn- is the characteristic function of 17^. If u G H{a,g) then 

[ii]r = a and Raniu) = g. 

Note that Hq{Q) can be identified with H{0, 0), and so the following defini- 
tion is consistent with the previous one. 

Definition 2.2. A function u G H{a,g) is a weak solution of ( ^.la| ), ( |2.1b ), 
(|2.1cD if V = u — g + axn- satisfies 

(2.6a) -B[v,^P]=F{i;), 

for all ip G Hq[^), where 

(2.6b) F(V) =Fi(V') + ... + F4(V') 

= y"/V' + JpVg-V^p- j [5Va-Vi) + j V ■ (b^jn). 
n n Q- n- 

A classical solution of ( |2.1a| ), ( |2.1bD , ( |2.1cD is necessarily a weak solution. 

Theorem 2.2. If f e L^{^), g G H '^j ^), a nd a, b e H^{^), then there 
exists a unique weak solution of ( \2. 6c\ ), ([2. 61^ , ( ^.Ic ) in H{a,g). 
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Proof. The right-hand side of ( |2.6b ) is a continuous hnear functional on 
Hq{Q). By the Riesz representation theorem, there exists a unique v £ 
Hf{n) such that -Blv,^;] = Fiip), for all ^ G ^^o(^)- Define the solution 
as u = V + g - axn- ■ □ 

Remark 1. We could replace / V • (bipn) by / bip ds in ( |2.6b| ). Then the 

n- r 

requirement for b is reduced to 6 G L-^(r), which is a bit weaker than our 
assumption ( ^^ ) for b. 

3. Abstract Form of Finite Difference Methods 

In this section an abstract framework is given for proving convergence 
of finite difference methods for the elliptic interface problem introduced in 
the previous section. 

For any h > 0, let Hq^ be a finite dimensional vector space with norm 
|[ • \\h. This space should be thought of as a discrete approximation of the 
Sobolev space H^i^) with grid size measured by the parameter h. 

On each finite dimensional space H^'^, we suppose there exists a bounded 
extension operator : Hq^ —>■ Hq{^}) with the bound 

(3.1a) \\T\i^'')\\Hi^^n)<CoU%, 

for all il)^ G ^o'^i with Cq independent of h. 

Furthermore, we assume a strong approximation property. That is, 
given ip € i/o(r2), there exists a sequence il^^ G H^'^ such that 

(3.1b) T''{^P^)^iP in H^{n), as h ^ 0. 

On each vector space Hq^, we suppose there is a uniformly bounded 
family of bilinear forms B^[-, ■] such that for every u^, G 

(3.2a) Ci\\u''\\l<B''[u^,u''] and \ B^[u'' ,v^] \< C2\\u^\\h\\v^\\h, 

for positive constants Ci, C2 independent of h. We also assume there exists 
a uniformly bounded family of linear functionals F^{-) on H^'^ such 
that for every ip^ G 

(3.2b) I F'^i^l^^) |< C3\\i;%, 

again with a constant C3 independent of h. 

Finally, we impose weak consistency with the weak problem ( ^.6a ), 
( 2.6b ). For every pair of sequences v^, G -^0''^ such that T^{v^) v 
weakly in H^i^) and T'^{ip'^) ip strongly in H^i^), we have that 

(3.3a) B^lv'',^;'']^ B[v,ij] 

and 

(3.3b) F'^(V;'»)^F(V^), 
where B and F are defined by ( |2.2| ) and ( p.6b ). 
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Under these general assumptions, we have the following: 

Lemma 3.1. For every h > 0, there exists a unique solution € H^'^ of 
the discrete problem 

(3.4) -B''[v'',i;''] = F^iij''), 

for every ijj^ E ^q'^. 

The sequence of extensions T^{v^) of the family of discrete solutions con- 
verges weakly to v = u — g + oxq- in Hq(Q), where u the solution of the 
weak problem ( |2.6a| ), ( |2.6b ). 



Proof. The existence of a unique solution of the discrete problem (3.4) fol- 
lows by the Lax-Milgram lemma. 

By the estimates (|3]2|), ( |3^ , we obtain a uniform bound for the se- 
quence of discrete solutions 

\\v%<C, 

with C independent of h. Thus, using ( p.la|) we have the bound 

\\T''{v')\\H^^in)<C. 

By weak compactness in the Hilbert space H^{n), there is a subsequence 
r'^ {v^ ) converging weakly to some v G Hq{^). 

Now let tp € Hq{Q) be given. Using the approximation property (3.1b), 
choose a sequence i/j'^ € H^'^ such that 

T^{il:^)^i) in Hl{VL). 

By weak consistency ( 3.3a ), ( |3.3b| ) we find that the limit function v sat- 
isfies d^ei), (|2!6b|) . Finally, since this problem has a unique solution, it 
follows that the full sequence T^{v^) converges weakly to v in Hq{Q). □ 



Remark 2. The simplest way in which to obtain a consistent, and there- 
fore converg ent, scheme is to define B'^[v^,^p'^] = B[T^{v^),T''('iIj'^)] and 
ph^^h^ _ i?(j'^(^^)). This is essentially the method of finite elements 
which then hinges on the choice of the extension operator T^, |||]. 

However, the scheme under consideration in the next sections is not of 
this type, insofar as it originates from different approximations for B and 
F. 



4. Numerical Method 

In this section, we rederive the finite difference scheme fro m [| | for ( ^.la ), 
( 2.1bD , ( |2.1c ), by discretizing the weak formulation ( 2.6a ), ( 2.6b ). 

For the remainder of the paper we assume that the coefficient matrix 
is of the form (3 Id, for some scalar function. We also assume that the 
data functions a, b, g, /, all lie in (7^(17), with a, b vanishing on dQ. For 
simplicity, we restrict ourselves to the special case of a rectangular domain 
O = {xw,xe) X {ys,yN) in the plane. Given positive integers / and J, 
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set Ax = {xe — xw)/{I + 1) and Ay = [un — ys)/{J + 1)) and define a 
uniform grid = {{xi,yj)} wliere Xi = xw + i^x and yj = ys + jAy 
for i = 0, 1, • ■ ■ , I + 1 and j = 0, 1, • • • , J + 1. The grid size is defined as 
h = min(Ax, Ay). The ratio Ax/ Ay is fixed when the grid size h goes to 
zero. 

The set of grid functions will be denoted by 
(4.1a) H^^'' = {w^ = (wij) : < i < / + 1, < j < J + 1}. 

The discrete solution space is defined as 

(4.1b) H^'^ = {ip^ = {ip.j) G H^^^ : ip^j = on the grid boundary}. 

To construct the bilinear form on Hq'^, we discretize the coefficient (3 in 
two ways as follows, 

(4.2) Pi+l/2,j = PiXi+l/2,yj), Plj+1/2 = Pi^h yj+1/2)- 

For E H^'^ define the usual finite difference operators 
(Va;w)i+i/2j = {wi+i,j - Wij)/Ax, 
yw)i j+i/2 = {wi,j+i - Wij)/Ay. 

For , G ^0^1 bilinear form then is given by 

J / 

(4.3) i?'^[^;^V^'^]=^J;/3^,/2^^.(V,^;)i+l/2,,■(V,V).+l/2JAxAy 

j=l i=0 

I J 

+ ^Y1 /5i+i/2(Vy'«^)ij+i/2(Vj;V')ij+i/2AxAy. 

i=l j=0 

As in the continuous case, we use this to define a norm on Hq'^: 

(4.4) \\i;^\\l = B^['ip'','ilj% 

Our next task will be to discretize the linear functional F in ( 2.6b| ). The 
data functions naturally give rise to grid functions a^, b^, etc., by restriction 
to the grid The normal vector is discretized by 

(4.5) n'^ = (n«'\n(2).'^) = (,/,^, </,^)/|V0'^| 
where for i = 1, • • • , /, j = 1, ■ • • , J, 

{(t>x)i,j = i(t>i+i,j - </'i-i,j)/(2Ax), 
i^y)i,j = - '^j,i-i)/(2Ay). 

Hence is defined on all interior grid points. It will not be used on 90. 
We can immediately define 

/ J 

(4.6a) F^ii,'^) = ^Y1 /^..j-^^j AxAy, 

i=i j=i 
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and 

(4.6b) F^i'il^'') = B%'',i;'']. 

The other two pieces require the locahzation of integrals to the subdomain 
First, we discretize the characteristic function: 



(4.6c) Xi,j 



1, if < 
0, if (pij > 0. 



Define 

(4.6d) xl+i/2,j = - 0i+i/2,j) + Xi+i,A+i/2,j) 

Q . ^ f + IfAi+ijI), if \4>i,j \ + \4>i+i,j\ > 0, 

*+i/2,i I Q otherwise. 

The factor Xi+i/2j^-^ approximates the length of the portion of the arm 
from {xi,yj) to (xj+i,yj) within 0". And also, let 

(4.6e) X-j+i/2 = - ^i,j+i/2) + Xi,j+i0ij+i/2) 

Q . ^ \\(t>i,j+l\l{\(t)i,j \ + l<^jj + l|), if \(t>i,j\ + \(t^i,'J + l\ > 0: 

«>i+i/2 I g otherwise. 

Now we define the remaining two pieces of the linear functional F. 
J I 

(4.6f) F^ii;'') =^Y1 V.(fcn(i)V)^+i/2,i AxAy 

j=l 1=0 

I J 

i=l j=0 

and 

J / 

(4.6g) F^ii;^) =-Y.Yl /?'+i/2,i(Vxa)i+i/2j(V,V)^+i/2j X'+i/2,, AxAy 

j=l j=0 

/ J 

/'?i+i/2(V2/a)ij+i/2(Vj,V')i,i+i/2 X-j+i/s^^Ay- 

i=l j=0 

Using (p^,. . . ,(|^), define 
(4.7) F'^ = Ff + 



With definitions ( |4.3| ) and ( |4.7| ) the discrete problem is formulated as in 

Next we show that the discrete weak formulation is the same as the finite 
difference scheme introduced in |^]. 
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Make the substitution = — + a^x^i using ( [4.6c| ), in ( ^ ) to write 

(4.8) - B%^, + B%^, - B\a^x^, i^^] 

We note right away that the second term on the left cancels with 
on the right. We are going to apply summation by parts to remove the 
difference operators from the test vector ip^. The idea is expressed by the 
one-dimensional formula 

/ I 

(4.9) - ^ ai+i/2(VxV')i+i/2 = xa)i'ilJi, 

1=0 i=l 

provided that = = 0. Here, {Vxa)i = (ai+1/2 - "j-i/2)/^3;- 
Using ([4.S|), the first term in ( [4.8D can be rewritten as 

/ J 

(4.10a) -B'^[n^ ^A'^] = 2^[V.(/?V.7z)ij + Vy{f5^Vyu)i,j]i^i,jAxAy. 

i=i j=i 

In the same way, we have using definitions ( 4.6d| ), ( |4.6e| ), 
/ J 

(4.10b) F^i^'') = -YYl + n'^jiVyX^Mj^xAy 

i=i j=i 

To treat the remaining term on the left-hand side of ( ^.81 ), we use the 
following product rule for the difference operator 

Vx(ax)i+l/2,i = {^xa)i+i/2jxl+l/2J + 0'i+l/2,jC^=^X)i+l/2,j 
'^i+l/2,j = - ^j+l/2,j) + 0.i,jGi+l/2,j-, 

in which Xi+1/2 j ^i+i/2,j were defined in (4.6d). Similarly, we have using 

(H) 

Vs/(ax)iJ+l/2 = Vyajj+i/2X?,j+l/2 + '^?j+l/2^a;Xij+l/2 
^'i,j+l/2 = ~ ^ij+1/2) + aj,i^j,j+l/2- 

It follows from this and ( [4.9D that 

i?'^[aV,V'"]+i^3'(V'") 
/ J 

= XlZl^^i/2,i«i+i/2j(V=«x)i+i/2,j(Vx^)i+i/2jAxAy 
/ J 

(4.10c) + E E /^L■+l/2«L•+l/2(VyX)^,,+l/2(V^,V)^,i+l/2AxAy 

i=l j=0 

/ J 
i=i 0=1 
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Combining (|l0|), dilObl) , (|4l0^ ), we obtain 

/ J 

i=i j=i 

I J 

/ J 

i=i j=i 
I J 

i=i j=i 

Since this must hold for all test vectors il^^ € Hq^, we have shown that the 
finite difference scheme 

holds at all interior grid points. Note that this is the scheme that was found 
in p, see equation (77) therein. 

Remark 3. Here the passage from the weak formulation via summation by 
parts to the finite difference scheme is analogous to what is often done with 
PDE's. 



5. Convergence 



In this section, we establish the converge of the scheme (3.4) 



Theorem 5.1. Let Q cM? be a rectangle. Assume that the data functions 
a, h, g, f all lie in C^{Q), with a, b vanishing on d^. Suppose that the 
coefficients have the form (3 Id. Then there is a family of linear extensions 

: Hq^ Hq{Q) which together with the bilinear forms (|4.3| ), and the 
linear functionals ( 47|) satisfy the structural conditions ( |3.1aD , ( |3.1b ), 
(|]2|), (|]2|), (|3|), ( pbl) . 

The sequence of extended approximate solutions T^{v^) of the discrete 
weak problem (|3.4D , ( [4. 3D , (|4.7j ) converge weakly to the weak solution of the 
PDE (|]6|), (ifebl ) in H^{n). 



Proof. The second statement follows from the first by Lemma 3.1. The next 



three subsections are devoted to the verification of the structural conditions: 
uniform boundedness, extension and approximation, and finally consistency. 

□ 
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5.1. Uniform boundedness. 

Lemma 5.1. The family of bilinear forms ( |4.3| ) satisfy the uniform bounds 

Proof. Because of our choice of norm in (|4.4|), the lower bound in ( ^.2a ) 
is immediate. The upper bound in ( 3.2a| ) follows easily from the Cauchy- 
Schwarz inequality. □ 

Lemma 5.2. The family of linear functionals ( [4.7| ) satisfy the uniform bound 
( pbD . 

Proof. We treat the four pieces F^, . . . , in turn. 
By the Cauchy-Schwarz inequality, we have 

The first factor is just a Riemann sum for the L^-norm of /, and so it is 
bounded by 2||/||2,2(f^), for h small enough. The second factor is estimated 
using the discrete version of the Poincare inequality 

/ J 

(5.1) (J^lV^^f A:EAy)V2<c(l^)(^^|(V,.^),+i/2,,f AxAy)i/2, 

i,j i=0 j=l 

which follows as in the continuous case using summation by parts. This last 
sum is then bounded by because of the uniform lower bound for 

the coefficients (3. 

Prom ( |4.6f ), we can use the discrete product formula to rewrite F^ as 
I J 

i=0 3=1 

I J 

1=1 3=0 

Since b and (n(^\ n^^)) are in C'^ and {x^ ^x^) is bounded, using the Cauchy- 
Schwarz inequality again, we get the bound 

\F'i{il^^)\ < C\\bnU2^n) U% + C||V • (HIIl^CH) | ^ M^^^^V ' 

Applying ( |5.1| ) to the last sum above, we obtain the desired bound for F^. 
In the same way, we get, 

1^2^^^^)! < C\\pVgh2(^^)\\^P% and \F^i^p'')\ < C\\(3Vah2^n)\\i^%. 

□ 
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5.2. Strong Approximation. In this section we define a uniformly bounded 
family of extension operators using the basic approach from the theory 
of finite elements. Then we verify the strong approximation property. 

To this end, for each h consider a triangulation of the domain J7 containing 
all triangles with vertices 

{{xi,yj), {xi+i,yj), {xi,yj+i)} or {{xi,yj), {xi^i,yj), {xi,yj^i)} 

based on the grid Q^ , as shown in Figure |5.2| . For any grid point (xj, yj) G 
ri'^, let rj^j E H^{Q) be the continuous, piecewise linear function which is 
equal to 1 at the grid point {xi,yj) and equal to at all other grid points. 
Given G H^'^, the extension operator is defined as 

id 

Then : H^'^ H^{n) and : h]'^ H^{^)- Although it is not 
required by the abstract framework, the extensions are linear operators. 

Since T^{w^) is the unique continuous linear interpolant on each triangle 
{{xi,yj), {xi±i,yj), {xi,yj±i)}, we have explicitly 

(5.2) 

T'^{w'^){x,y) = Wij + {Vxw)i±i/2j{x - Xi±i) + {Vyw)ij±i/2{y - yj±i). 
Given the coefficient function /? on J7, we write using the definition ( [4.2| ) 

In general, given a discrete function 13^'^ defined at the half grid points 
(xj_|_i/2) yj)) s-s above, we define the piecewise constant extension 

S''\p'''^){x,y) = Pl,/,^^ 

on every triangular cell having the horizontal edge from (xi,yj) to (xj+i, y^), 
see Figure 2. Similarly, given a discrete function /3^''* defined at the half grid 
points (a^j, yj+1/2)) we define 

S'''^{P'''^){x,y) = Pl^^y, 
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Figure 2. 



on every triangular cell having the vertical edge from {xi,yj) to {xi,yj+i), 
see Figure 2. 

The key to our estimates will be to replace sums by integrals. The fol- 
lowing lemma summarizes the important formulas. 

Lemma 5.3. With the definitions above, we have 



(5.3) B^[v^,iIj^] 



(5.4) F^ii;^) 



(5.5) F^ii;^) 



J 

Jn 



+ S^'''{P'^'^)S'^'^{x^'^)T^{a'')yT''{ij;'')y 



Proof. Suppose that j '^^ defined for i = 0, . . . , /, j = 0, . . . , J + 1, 

and ^l_^_i/2 J = for j = or j = J + 1. Each point of the form (xj+1/2, Vj), 
i = 1, . . . ,1, j = I, . . . , J, lies on the boundary of the two triangular cells 
sharing the arm from {xi,yj) to {xi+i,yj) , and on these cells we have 
^i+i/2j ~ ^^'^{C^'^){^iy)- On the remaining cells, namely those with j = 
OT j = J + 1, we have ^l^i/2j ~ S'^ 'H'^^ '') = 0- Therefore, since the area of 
each pair of cells (where S^'^{^^'^) could be nonzero) is equal to Ax Ay, we 
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have 



i=0 j=l 

Similarly, we have 

/ J 



i=l j=0 

provided that = 0forz = 0ori = I + l. 

Notice that the piecewise continuous extensions S°''^ are multiplicative 
in the sense that, for a = 1,2, 5°''^(a"''*6"''') = 5"'''(a"''^)5"''*(6"''*), for 
arbitrary a°''^, h°^'^ defined at the half grid points. So using the formulas 
just derived, we have from (^]^) that 

Jn 

+ S^'''{^^''')S^'^{VyV'')S^'^{Vyi^'')]. 

By Q, we see that S^'^{V^^p^) = T^{il)^)x and S'^^^{Vy^^) = T^{'^^)y, 
and therefore, we have verified ( |5.3| ). 

The proofs of the other two formulas are similar and will be omitted. □ 

We are now ready to study the extensions. 

Lemma 5.4. The family of extensions : Hq^ — > Hq{Q) satisfy the uni- 
form estimates 

j^Wi' \\h < U V' < \\h 

Proof. Let il^^ E ^o'^- By definition, we have 

||rVllHi(f7) = B[T^i,^ ,T^ij^] and Wi^^Wl = V'']- 
Moreover, by (|5.3|), we have that 



Now by the uniform bounds on /3, we have that 

(3{x,y)<—(3{x',y'), 
m 

for arbitrary {x,y), {x',y') G fi, and so 

(5.6) < S^'^iP'^^), S^'^iP^'^) < —p. 

M m 
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By (|]D, it follows that 

m M 
M m 

□ 

Lemma 5.5. The extension operator : Hq'^ Hq{Q) has the strong 
approximation property ( 3.1b| ). 

Proof. First let tp £ Cq{^1), and define ijj^ G Hq'^ to be the grid function 
with values V'jj' = i^{xi,yj)- Since ip G C^{Q), ip and its first derivatives are 
uniformly continuous on 17. It follows from ( |5.2| ) the piecewise constant func- 
tions T^{tP^)x and T^{il)^)y converge uniformly to il^x and tpy, respectively, 
therefore also in L'^{uj): 

||V(rV-^)llL^(n)-0, 
as h ^ 0. Since P has a uniform lower bound, this implies that 

as — s- 0. 

Finally, the result for general tp G -frg(r2) follows by density. □ 
5.3. Weak consistency. 

Lemma 5.6. The bilinear form ( [4.3D satisfies the weak consistency hypoth- 



esis ( 3.3a ). 

Proof. Suppose that v^, ip^ G Hq'^ are sequences such that T^{y^) v 
weakly in Hq{Q) and T^{'il)^) —>■ ip strongly in H^{Q,). Thus, we have that 
VT^{v^) Vv weakly in L'^{n) and VT\ip^) Vip strongly in L'^{n). 

Recalling the definitions from the previous subsection, we have that the 
functions {a = 1,2) are both uniformly bounded and converge 

pointwise to /3 in \ F. Writing 

S''\p'''')T^{ij% - Pi^x = [5i''^(/3i''^) - /3]V. + S^'^{(5^'^)[T^['P% - V^,], 

it follows by the dominated convergence theorem that S^^^{[3^^^)T^{'4)^)x ^ 
fiipx strongly in L^(Q). 

Now, using Lemma ( |5.3| ) we have that 



B'^[v^,^'']= I [S^'''{(3^''')T''{v^)xT''{iP^) 

+ S'^^^{0^^^)T^{v\T^{'4j^ 

'^l.h ( nl,h\'-rh f „i.h 



= {T\v%,S''\(5'^'')T\i^^)x)LHn) 

+ {T\v%,S'^\(3'''^)T\^%)l.^^) 

= B[v,i>], 

as h^O. □ 
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Lemma 5.7. The linear functional F defined in (O) satisfies the weak 
consistency hypothesis ( |3.3bD . 

Proof. Let € -f^o(^) aiid suppose that G Hq'^ is a sequence such that 
T^{il;^) V strongly in Fd(J^). We must show that F^{ip^) F{'il)). 

To begin, we observe that it is enough to prove this for test functions in 
Cq{Q). For any € Cg(r2), we can write 



\F^{'4,^)-F{'4,)\ < iF^i^p'' 

< iif'^iiiiv'^ 

By Lemma ^.4|, we have 



fy + \FHA 



l.h 



m 



< 



M 



m 



■i^'')\\H^{n) 



By Lemma 5^, the norms 
we have shown that 



are uniformly bounded. Therefore, 



We claim that this can be made arbitrarily small for all h < h^. By as- 
sumption, we have that T^{'>p'^) in Hq{Q,). By density, tp can 
be chosen arbitrarily close to tf^ in Hq{Q). By construction, we have that 

T^{ip'^) — -0 ^ in Hq{Q). This covers all but the first term above. The 
first term can also be made small if the consistency condition is valid for 

We now proceed to verify the consistency of F under the assumption that 
ip G Cq{Q) by considering each individual piece. 

By ( [4.6a ) and the fact that ip S Cg(r2), we see that F^^iij:^) is simply a 
Riemann sum for Fi(V'), and thus (0^) ^ Fi(V'), as /i ^ 0. 

Given g G (7^(0), we have that T^{g^) ^ g in Hq{Q), as in the proof of 



Lemma Thus, we may apply the result of Lemma p.q to conclude that 
F^ii^^) = B%\ij''] ^ B[g,ij] = F4(V). 



According to (5.4), we have 



[^i,h(^i,h);^h(^h^{i),h^h)^ + S^'''{x^''')T''{b''n^^^'''ip\]. 



Since 6, n, ^/^ G (7-^(0), we have as in the proof of Lemma 5.5, that VT^(6^ra^'0^) 
Vbnip in L^. Moreover, 5°''^(x"''^) is uniformly bounded and tends to Xn- 
pointwise in Q \ F. Thus, it follows that F^{il:'^) F^tp). 

The consistency of F^ follows from the formula ( |5.5| ) , in the same way. □ 
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